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Abstract 

The pionium decay width is evaluated in the framework of chiral perturbation theory 
and the relativistic bound state formalism of constraint theory. Corrections of order 0(a) 
are calculated with respect to the conventional lowest-order formula, in which the strong 
interaction amplitude has been calculated to two-loop order with charged pion masses. 
Strong interaction corrections from second-order perturbation theory of the bound state 
wave equation are found to be of the order of 0.4%. Electromagnetic radiative corrections, 
due to piou-photon interactions, are estimated to be of the order of —0.1%. Electromag¬ 
netic mass shift insertions in internal propagators produce a correction of the order of 
0.3%. The correction due to the passage from the strong interaction scattering amplitude 
evaluated with the mass parameter fixed at the charged pion mass to the amplitude eval¬ 
uated with the mass parameter fixed at the neutral pion mass is found to be of the order 
of 6.4%. 
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1 Introduction 


The 7T7T scttering amplitude [ill] represents one of the main quantities that allow confronta¬ 
tion of predictions of chiral perturbation theory 031 with experiment. Unfortunately, 
the lack of direct low energy data forces one to reconstruct the low energy scattering 
amplitude from extrapolations P, ^j. |7j from high energy data |j. |l^] and from indirect 
information coming from decay 0, at the price of increasing error bars on numeri¬ 
cal values. The 20% uncertainty of the “experimental” value of the isospin zero S'-wave 
scattering length, a® = 0.26 ± 0.05 1 0 - 0 - 0 , does not allow one to draw a clear-cut 
conclusion when the latter is compared with the theoretical prediction of standard chiral 
perturbation theory, which is 0.20 to the one-loop order || and 0.217 to the two-loop 
order 0- 

From this viewpoint, the DIRAC experiment, which will be realized at CERN in the 
near future and which aims at measuring the pion scattering lengths from the lifetime of 
pionium (7r + 7r~ atom) decaying into 7r°7r° ifld] . 0, might provide a decisive improvement 
for the above comparison. (Results on previous experiments are presented in Ref. Hl8| .) 

On the other hand, it was noticed |0 on theoretical grounds that the standard for¬ 
mulation of chiral perturbation theory does not provide the most general basis for an 
interpretation of the experimental results within the framework of QCD. It is generally 
assumed that the quark condensate in the chiral limit, — < 0|gg|0 > 0 /iy?, has a mass 
scale of the order of 1 GeV, typical of massive hadron masses [[5|. 0. This hypothesis, 
while apparently natural, has not yet received a direct experimental confirmation. In this 
connection, a generalized form of chiral perturbation theory was proposed 0- 0 


m 


which the above hypothesis is relaxed and the quark condensate itself appears as an ad¬ 
ditional expansion parameter, allowing, in certain processes, its experimental evaluation. 
Thus, the tht scattering amplitude A(s\t,u) becomes at lowest order 0(p 2 ): 


A(s\t,u) = —^(s-2rhB 0 ), 


Bo — — < 


>o /Fl m = 


m u + m d 


( 1 . 1 ) 


F 0 is the pion decay constant F n in the chiral limit (F n = 92.4 MeV), m u and m d are 
the masses of the quarks u and d; Zweig rule violating effects have been neglected. The 
quark condensate parameter 2 7nB 0 is expected to take values between 0. and , the 
latter value reproducing the standard predictions of chiral perturbation theory. 
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For 2fnB() ~ 2 m i f° r instance, one finds for the scattering length a[ j, including one- 
loop effects |^2], |23|], the value 0.27, which lies 35% away from the standard value 0.20 and 
closer to the experimental value 0.26. This example underlines more acutely the necessity 
of disposing of precise experimental informations on the pion scattering lengths to be able 
to distinguish between various chiral symmetry breaking schemes. 

Coming back to the pionium lifetime, its lowest-order expression was established long 


ago in the nonrelativistic limit by various methods [24]. £5], [26]. [27 
1 167T /2A m n (a° — Oq) 2 


— = T 0 = 

To 


9 


m-n 


m 


'|V ? +-(0)| , Am n = m n + — m^o, 


( 1 . 2 ) 


where a J 0 is the (dimensionless) S'-wave scattering length in the isospin / channel, usually 
evaluated in the literature with the charged pion mass, and ^.(O) is the wave function of 
the pionium at the origin (in axspace). Characteristics of the pionium have been discussed 
in Ref. ^8| and the relevance of its lifetime for determining chiral symmetry breaking 
parameters has been outlined in Ref. ]^9|. While the above formula provides a relationship 
between the pionium lifetime and the pion scattering lengths, it is desirable, for a precise 
theoretical interpretation of the experimental result, to have a knowledge of the possible 
corrections to it. This question was addressed recently in Refs. [30|. |3T], [32]. |33| . 

The purpose of the present paper is to evaluate the corrections to formula ( |1.2|) in the 
framework of SU (2) x SU (2) chiral perturbation theory. Apart from relativistic kinematic 
and mass shift corrections, they can be grouped into four categories: i) Corrections coming 
from second-order perturbation theory in the bound state wave equation, ii) Contributions 
originating from the electromagnetic radiative corrections due to pion-photon interactions, 
iii) Contributions coming from the electromagnetic mass shift corrections, due to quark- 
photon interactions and acting through insertions of the 0(e 2 p°) mass shift lagrangian 
term in pion internal propagators, iv) Mass shift corrections with respect to the strong 
interaction amplitude evaluated with the charged pion mass. 

The evaluation of the pionium bound state energy shift is done in the framework of 
constraint theory relativistic wave equations |34], which can be considered as a variant 
of the quasipotential approach [35|, [36] and has been shown to provide a means of a 
covariant treatment of the QED bound state problem [37|. The above corrections to the 
bound state energy shift are evaluated to the relative leading order 0(a), where a is the 
QED hne structure constant, the calculations being generally done to one loop (globally 
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with respect to the strong and electromagnetic interactions). 

Our results are the following. The corrections of the first type are found to be of 
the order of 0.4%. The corrections of the second type are shown to be free of infra-red 
enhancement and are estimated to be of the order of —0.1%. The corrections of the third 
type are estimated to be of the order of 0.3%. The corrections of the fourth type are 
found to be of the order of 6.4%. 

The plan of the paper is the following. In Sec. 2, the properties of the constraint 
theory bound state equation are briefly sketched. In Sec. 3, the latter formalism is 
adapted to the case of the coupled channels of the 7r + 7r - and 7r 0 7r° systems. In Sec. 4, 
the pionium lifetime expression in first-order of perturbation theory with respect to the 
strong interactions with pion mass shift is established. In Sec. 5, the corrections due 
to second-order perturbation theory in the bound state equation are evaluated. In Sec. 
6, the radiative corrections due to the pion-photon interaction are evaluated to one-loop 
order in the tree appriximation of the strong interactions. In Sec. 7, the electromagnetic 
mass shift corrections are calculated. A summary of the results is presented in Sec. 8. 
Several details of calculations are presented in Appendices A-E. 


2 The constraint theory bound state equation 


The Bethe-Salpeter equation |38|], which is the basic bound state equation in quantum 
held theory, has been revealed inadequate for quantitative calculations with covariant 
propagators. Two typical drawbacks are the following. In the nonrelativistic limit, the 
one-photon (or one-particle) exchange diagram yields relativistic corrections of order 1/c, 
instead of 1/c 2 |IU|. In spectroscopic calculations, two-photon exchange diagrams yield 
spurious infra-red logarithmic singularities J40j] . These effects are cancelled only with the 
inclusion of higher order diagrams, a feature that enormously complicates the use of the 
equation in perturbation theory. 

In practice, the Bethe-Salpeter equation has been used in QED in the Coulomb gauge, 
which is a noncovariant gauge. Because of the instantaneous nature of the dominant part 
of the photon propagator, one is able to transform the original four-dimensional equation 
into a three-dimensional one and to avoid the previous difficulties ||I[ . However, the latter 
gauge has its own limitations. It neccesitates a different treatment of exchanged photons 
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and of photons entering in radiative corrections. Furthermore, additional complications 
arise when QED is mixed with other interactions, where already covariant propagators 
are present. 

In this respect, the wave equations obtained in the framework of constraint theory 
|34j have been shown to provide a satisfactory answer to the requirement of a covariant 
treatment of perturbation theory in the bound state problem |5j] . 

In constraint theory the relative energy variable is eliminated by means of a constraint 
equation. For a two-particle system this is generally chosen in the form: 


C(P.P) = (p, - p 2 ) - (m J - m 2 ) fs 0. 


( 2 . 1 ) 


The two-particle Green’s function is projected on this hypersurface and then iterated 
around it. At its pole positions, the latter not being affected by the projection operation 


(|2 . 1|) , one establishes a three-dimensional eigenvalue equation that takes the form: 


g~H = -VV, 


( 2 , 2 ) 


where gf 1 is the wave equation operator, which will be specified below, and V is the po¬ 
tential, related to the renormalized off-mass shell scattering amplitude T by a Lippmann- 
Schwinger type equation: 

V = T + Vg 0 T, f = ^=T\ c , (2.3) 

where the index C denotes the use of constraint ( |2.1| ) (on the external lines of T) and s = 
(Pi+Pi) 2 - (T is defined as the amputated four-point connected Green’s function multiplied 
by the wave function renormalization factors of the external particles.) The amplitude T 
contains the usual Feynman diagrams, where the external particles are submitted to the 
constraint C. The second term in the right-hand side of the first of Eqs. (|2.3|) generates 
an iteration series, the diagrams of which are called “constraint diagrams”, where the 
integrations, because of the presence of the factor g 0l are three-dimensional, taking into 
account constraint C. 

As long as perturbation theory is concerned, Eq. (|jjj) is equivalent in content to the 
exact Bethe-Salpeter equation, with, however, a different arrangement of the perturbation 
series. The constraint theory wave function T is related to the Bethe-Salpeter wave 
function $ by means of the projection of the latter on the constraint hypersurface (|2.1| ), 
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but the explicit form of this relationship will not be needed in the present work. Once 
Eq. ( |2.2| ) is solved with the exact potential V, the Bethe-Salpeter wave function $ 
can be reconstructed, through the iteration procedure, in terms of the constraint theory 
wave function T and can be shown to satisfy the Bethe-Salpeter equation with the exact 
irreducible kernel K and the same energy eigenvalue as that of Eq. ( p.2|) |37 | . 


The choice of the operator g 0 is not unique in principle, but in practice it is made 
on the basis of several natural requirements. The choice that is made below satisfies the 
following four properties: i) Correct nonrelativistic limit (Schrodinger equation), ii) Cor¬ 
rect one-body limit. When one of the masses becomes infinite, one recovers the Dirac or 
Klcin-Gordon equation in the presence of the static potential, iii) Correct hcrmiticity and 
unitarity properties. Potential V is an irreducible kernel, in the sense that it is free of sin¬ 
gularities in the s-channel, at least in the elastic unitarity region: the constraint diagrams 
cancel the singularities of the reducible diagrams of T. iv) Correct QED spectroscopy. In 
particular, the constraint diagram contributions remove all spurious singularities (in the 
bound state region) coming from T. More generally, the leading effect of the sum of all 
n-photon exchange diagrams (in the absence of radiative corrections) is of order 0(ct 2n ), 
where a is the fine structure constant. 

When constraint ( [2.1D is used, the Klcin-Gordon operators of particles 1 and 2 become 
equal: 


H 0 = {p\ - ml)\ c = (p| -ml) | c = m\ + m 2 ) + — 

We use the notations: 


4P2 


m 2> +J) T2 


(2.4) 


P = Pl+P2, P=^(P1~P2), X = i(xi + X 2 ), X = X\ 


X‘2 j 


(2.5) 


and the decompositions of four-vectors into transverse and longitudinal vectors with re¬ 
spect to P: 


<ll = % 


pP p 

p2 r 


<£. = (q-p)p^ 


p _ m 


ql = q-P , 


Pl = v/P 2 , 


r — \ —x 


r T 2 


( 2 , 6 ) 


For two spinless particles, go is chosen to be: 

1 

9o = 


Hn + ie 


(2.7) 
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up to possible finite renormalizations related to the finite parts of the individual parti¬ 
cle propagator renormalizations; they will not show up, however, to the approximations 
used throughout this work. Furthermore, go undergoes a finite multiplicative renormal¬ 
ization by a constant factor (1 + 71 ) due to the off-mass shell treatment of the Lippmann- 
Schwinger type equation (|2.3| ). The constant 71 appears from the requirement that the 
only 0(1 jr) terms in the QED potential come from the one-photon exchange diagram. 
Its presence amounts to multiplying the potential V by (1 + 71 ) and continuing the use 
of expression (|2.7|) for g 0 , the constant 71 allowing the cancellation of a spurious 0(a 3 ) 


term j37|. The use of this finite multiplicative constant, which tends to improve the per¬ 
turbative expansion of the potential, should not, however, have an influence on physical 
quantities (in analogy with the presence of wave function renormalization constants). 

Since in the present work we are interested by corrections of order 0(a) to the pio- 
nium bound state energy, we can from the start consider the pure QED potential in its 
nonrelativistic limit (Coulomb potential) and use the corresponding nonrelativistic wave 
functions for the zeroth-order approximations. The pure QED corrections in the channel 
7r + 7r - — 7r + 7r _ , being of order 0(a 2 ) |37], [42|, will not be considered further. 

The nonrelativistic Coulomb potential is here: 


a 

Vcoui. 2 n , 
r 


h = 


17111712 

mi + m2 


( 2 . 8 ) 


The rest of the potential in Eq. (T3) will be treated as a perturbation. It contains 
the strong interaction part of the tttt interaction, as well as the interference part between 
strong and electromagnetic interactions. 

I 11 general, potential V being energy dependent, the scalar product of wave functions 
has a more complicated kernel than in the energy independent case |37| . The perturbation 
theory formulation in the case of energy dependent potentials can be found in Ref. |3(| 
(valid for four- and three-dimensional equations). However, since in the present work 
the zeroth-order potential is the energy independent Coulomb potential ( [2.8|) , the scalar 
product that should be used in the perturbative calculations is the usual nonrelativistic 
one. Energy factors, present in higher-order potentials, should then be expanded around 
their zeroth-order values. 

In the rest of this work we shall use, for the evaluation of the importance of various 
terms, the infra-red counting rules of the QED bound state system. Let, for a given 
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process 1 + 2 —> 3 + 4, s, t , m be the Mandelstam variables: s = (pi +P2) 2 , t = (pi — P3) 2 , 
u = (pi — Pi) 2 . We also define the (c.m.) momentum operators: 


b 2 ab ( 


T 2 l \ S 1/2 I 2 \ , 

K ,06 (s) = T--(m ff a+m fffc ) + — 


„T 2 


= -P 

-^) 2 


(a, b = 1, 2 or 3,4), 


4s 


= pIl ~ m l« = PbL ~ m l>‘- ( 2 - 9 ) 


In the pionium state (7r + 7 r _ ), the deviation of s from the threshold value 4m^+ is of order 
0(a 2 ). The quantities b+_(s), 6g 0 (s), 6 g + _(s), t and u are of order 0(a 2 ). The quantity 
^ 0,00 °f l^ ie 71-071-0 system at the same energy is of order 0{Am n /m. K ). 
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Wave equations of the 7r + 7r and 7r°7r° systems 


In order to deal with the specific sectors of the 7r + 7r“ and 7r°7r° systems, we have to 
enlarge the spaces of potentials and wave functions considered in Sec. 2. We introduce a 
two-component wave function T as: 





(3.1) 


and define the potential V in matrix form in the corresponding space: 

/ 


V = 


y boo,+- Km 


K—,+- K__ 0 o 

/ 00,00 

The constraint propagator g 0 [Eqs. (|2.2|) , ( |2.3|) and (|2.7| )1 is now: 

/ 


(3.2) 


\9o] = 


9o,+- 0 

0 1 


[dor 1 = 


( —1 

%,+- 


Po,oo 


0 


y 0 2p 0 


— 1 
00 


(3.3) 


where the subscripts H— and 00 have been associated with g 0l in order to take account, 
when necessary, of the specific free masses of the systems 7r + 7r^ and 7r°7r° in formula (|2.7| ). 
The factor 1/2 in front of g 0 ,oo l ias been introduced because of the identity of the particles 
in the corresponding sector. 

Let T + _ + _, T + _ j0 o, T’ Ci o :H _and T 00 00 (with M. = — iT ) designate the scattering ampli¬ 

tudes of the processes 7T + 7r~ —»■ 7r + 7r _ , 7r°7r° —> n + 7i~, etc.. We use for our calculations the 
chiral effective lagrangian P [|, |2l| in the SU(2 ) x SU(2 ) case. The scattering amplitude 










is obtained from the amputated four-point Green’s function of pseudoscalar densities, 
multiplied by the corresponding wave function renormalization factors. We shall use for 
the held U of the chiral effective lagrangian the representation 


7r 


U = a + i—.T, 

F 0 


7T- 




p 2 ' 


(3.4) 


where r are the Pauli matrices and 7r the pion helds. The pseudoscalar densities are 
defined as P a = iq^r a q (a = 1,2,3), where q are the quark helds. Their Green’s func¬ 
tions are obtained by deriving the generating functional with respect to the pseudoscalar 
sources. In the representation ( |3.4| ) and for the 0(p 2 ) lagrangian the pseudoscalar den¬ 
sities P a are proportional to the helds 7r a by the common constant factor 2B 0 F 0 . In 
the higher-order lagrangian terms more complicated differences arise and they should be 
taken into account. We shall draw Feynman diagrams with respect to the pion helds ix a . 

In the strong interaction limit and in the absence of isospin breaking, the amplitudes 
T above are related to the conventional amplitudes A(s\t,u), etc., with the relations: 


■rpstr. 

L1 +-,+- 

= A 4 + i+ _ = \t, u) + A(t s, u), 


o 

o 

-b r 

• <s> 

1 

= ^T^ + _ = Mf- fi0 = A(s\t,u), 


■rpstr. 

11 00,00 

= -^00,00 = u ) + u ) + A(u\t, s ), 

(3.5) 


and in terms of the isospin invariant amplitudes T iT \ they are 

T str. 

j_ 


= Ijd 2 ) i Ijd 1 ) _i ijdo) 
6 2 3 


pstr. 

-+-,oo 


rristr. 

1 00,00 


_ rpstr. _ 

°° !+ _ - 3 3 

_ ^ jd 2 ) i Ijdo) 

3 3 


The scattering lengths are defined as: 

— iT^^s = 4 ml) = 32naj =0 . 


(3.6) 


(3.7) 


Equations (|2.3|) take now the explicit forms: 


V+-,+- 

— + V+-,+-9o,+-T+-,+~ + -I+_ i oo5 , o,oo+o,+- 5 


J 

o 

o 

J 

O 

o 

+ 

+ 

1 

+ 

1 

o 

o 

+ 

to | I—‘ 

J 

o 

o 

© 

o 

o 

o 

o 


hoo,+- 

— T 00 , + _ + Voo,+-fi , o,+ 2 boo, 00 ^ 0,00 Too,+-) 


boo,00 

O 

o 

o 

+ 

o 

+ 

+ 

1 

1 

o 

o 

+ 

tO | H- 1 

o 

o 

o 

o 

o 

o 

o 

o 

(3.8) 


9 



We isolate from the potential f+_ + _ the Coulomb potential part [Eq. (|2.8|)] : 


Vh —— — Vcoul. + V 4—,+-• 


The wave equations (|2.2| ) then become: 


— {Vcoui. + ^+-,+-)^+- + f / +_ ) oo'I / oo, 


i x x 

— 9o, oo^oo = ^f / oo,+- v I / +- + -Eoo.oo'J'oo- 


(3.9) 

(3.10) 

(3.11) 


These two wave equations are characterized by the same eigenvalue P 2 . Whenever not 
specified, all potentials will be considered at the pionium ground state c.m. energy. At 
zeroth order (Coulomb potential only), it is: 

771/7T+ O 


Po = 2m, 


-a 


(3.12) 


For this value of P 0 , the operator g 0 ,oo hi Eq. (|3.11|) has values in the scattering region 
of the 7r°7r° system and therefore it gives rise to a scattering state for the wave function 
Too. The energy shift of the pionium is obtained by first eliminating Too from Eq. ( |3.10| ) 
in terms of T + _ through Eq. ( |3.11| ), using there the boundary condition that T 0 o is an 
outgoing wave due entirely to the presence of T + _. 

The wave function T 0 o can be expressed in terms of T + _ by iterating Eq. ( |3. 1 ID 
with its last term (proportional to Voo,oo)- This iteration, where the dominant contri¬ 
bution comes from the strong interaction sector, when treated globally as a first-order 


perturbation, yields in the pionium decay width expression the unitarity factor [^5| |27 
(1 + (2/9)(Am 7r /m 7r )(2aQ TOq) 2 ) -1 ; the correction term to one, being of the order o: 
will be neglected in the following. We thus obtain the wave equation: 


-4 


%X- xI> +- = 


Vcoul. + y+-,+- “ ^+-,00,90,00^00,+- 


T + _. 

(3.13) 

V*-,+- 

V +-.00 and 


Voo,+- are, in x-space, three-dimensional delta-functions and hence they project all mul¬ 
tiplicative quantities on their values at the origin. The value at the origin of the function 
g 0 ,oo is calculated by Fourier transformation to momentum space and dimensional regular¬ 
ization. Designating by A m n the pion mass difference [Eq. ( jl ■ 2| )j , one finds (see Appendix 
A, Eq. Q): 

(3.14) 


30,00 {r = 0) = — ^ Am n (m n + +?«#). 
This term induces an imaginary part to the pionium energy. 
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4 Lowest-order formula with mass shift 


We shall determine a lowest-order expression for the pionium decay width with the in¬ 
clusion in it of the main part of the pion mass shift. This formula will prove useful for 
the evaluation of the various types of correction. It is obtained by treating the last term 
of Eq. (|3.13|) in first-order of perturbation theory and by keeping in the potential Voo,+_ 
its dominant part, which comes essentially from the strong interactions. The expression 
of Voo,+- [Eq. (|3.8|) 1 contains three terms. The second and third terms are the constraint 
diagram contributions, which will be considered below. The Erst term is the scattering 
amplitude. The latter can be split into two parts, the strong interaction part and the 
rest, which represents the contributions containing electromagnetic effects. Among the 
latter, there is one piece which plays a crucial role; it is the quark-photon (or massive 
hadron-photon) interaction term in the chiral limit, which is reponsible of the main part 
of the pion mass difference |T|. It corresponds to the 0(e 2 p°) term of the chiral effective 
lagrangian with the expression e 2 C < QUQU 1 >, where C is an unknown constant and 
Q is the quark charge matrix [Til. [45| . Using for the held U the representation (|3.4|) , one 
Ends: 


2 C < QUQU ] >= 


2 e 2 C 


-7T 7T 


(4.1) 


which shows that this term induces for the charged pions the mass shift: 


= 2e 2 4p (4.2) 

which is nonvanishing in the chiral limit. Otherwise, it has no eEect on the scattering 
amplitude in lowest order. (We emphasize that expression O of the latter is also 
valid oh the mass-shell and does not depend of any mass-shell prescription.) Therefore, 
the quark-photon interaction term ( fl.l| ) acts essentially through insertions in the pion 
loop propagators, where the charged pion masses are replaced by their (almost) physical 
masses. Because the pionium lifetime evaluation is sensitive to the pion mass cliflerence, 
it is natural to incorporate from the start the quark-photon interaction term ( |4.1|) and its 
counterterms of the higher-order lagrangian in the strong interaction lagrangian. With 
this prescription, the amplitude A4 is split into two terms: 


M = M str - +qi + M em ', 


(4.3) 
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where J\A em ■ contains all interference terms between electromagnetic and strong interac¬ 
tion effects, except the lowest-order quark-photon interaction term and its counterterms, 
which are included, together with the strong interaction terms, in TVd str ' +<?7 . 

We now turn to the evaluation of the two constraint diagrams corresponding to the 
last two terms of the expression of the potential Voo,+- [Eq. (|3.8| )1 . The evaluation of 
these diagrams is done by first considering the potentials and amplitudes T str ■ at the tree 
level (Fig. 1) and calculating the loop integrals with the physical pion masses (cf. the 
comments above). The details of the calculations are presented in Appendix A. It is found 
that the terms in the amplitudes proportional to t, u , b 2 and 5 q yield, after integration, 
terms that are again of the same order and hence can be neglected. The leading terms 
are directly obtained by ignoring the above terms in the amplitudes and integrating over 
g 0 . One finds: 

T + -,+-g 0 ,+-T+-, + - = 

T + _,oo£o,ooToo, + - = 

(there are no integrations in the right-hand side,) and so forth for the other amplitudes, it 
being clear that go +- yields a real contribution and ( 70,00 an imaginary one. By the very 
choice of the constraint propagators go, the imaginary terms that arise from the constraint 
diagrams cancel similar terms coming from loop diagrams contained in the amplitudes T 
(the first terms of the right-hand sides of Eqs. (|3.8|) ). Therefore, these amplitudes are 
replaced by their real parts. Concerning the real part of the constraint diagram, it is of 
order 0(a ) and has the opposite value of the 0(a) part of the deviation of the scattering 
amplitude from the 7r + 7r“ threshold down to the pionium energy; it comes from the finite 
non-analytic part of the unitarity loop, the polynomial parts giving only contributions of 
order 0(a 2 ) to the above deviation. Hence, the real part of the constraint diagram shifts 
the real part of the scattering amplitude to its value at the 7r + 7r“ threshold. We thus 
obtain: 


y(°) 

v +-,+- 

= Relf2, + _(s = 4m^), 


y(°) 

> / oo,+- 

= Kefix? + -(s = 4mJ + ), 


y(°) 

v 00,00 

= Kef£] m (s = 4 mj + ). 

(4.5) 


The above property can also be generalized to the two-loop level of the strong inter- 


~^\/- 6 o,+-(r+-,+-) 2 , 

+ ^ \! ^ 0,00 ^+-,00^00,4— 5 (4-4) 
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action amplitude. The details of the derivation are presented in Appendix B. 

The part of the potential Voo,+- that will contribute to the lowest-order expression of 
the decay width is then provided by the real part of Adoo,+- 7 through TZeT^ + _. Defining 

T 


/h = /( 


OR 


' 2 ’ 


(4.6) 


one finds the modified lowest-order expression of the pionium decay width: 




Am, 

2r?VH 


(4.7) 


where IZeM^ '^V is calculated at the 7r + 7r threshold. (0 is the relative motion part of 
T.) Furthermore, because the external particle momenta are subjected to the constraints 
Pil = P2L = Pl/2, Psl = Pal = Pl/ 2, with Pi = 2 m n + and the kinetic energy operators 
\p 72 \ and \p' T2 \ have, in the bound state, values of the order of magnitude of 0(a 2 ), the 
mass-shell conditions in the above amplitude are: 


22222 a 2 i n 

Pi =P 2 = Ps =Pa = ? 7 br+; s = 4m^+, t = U = 0. 


(4.8) 


With respect to the lowest-order formula ( |1.2|) , formula (T7) contains two corrections. 
The first one, which has a kinematic origin, is included in the square-root term. The 
second one is included in Ad str +97 , where now, because of the pion mass difference, a 
modification occurs from the expression involving the scattering lengths calculated with 
the strong interaction amplitude. 

In the remaining part of this paper we shall evaluate the 0(a)-corrections to the 
formula (|4.7| ) as well as the modifications contained in _M. str - + i'y. 


5 Second-order perturbation theory 

In this section, we evaluate the effects coming from the second-order perturbation theory 
treatment of the strong interaction potential. At this order, it is the interference of the last 
two potential terms of Eq. (|3.13|) that contributes to the imaginary part of the energy. We 
have to distinguish here between the contributions of the discrete and continuum states 
of the pionium spectrum. We first consider the discrete spectrum contribution. 


13 











5.1 Contribution of the discrete spectrum 


Designating by "0+-,n (n > 1) the zeroth order radial exitation wave functions of the 
ground state [ij) is the relative motion part of 4/] and by E n the corresponding nonrela- 
tivistic energies, the shift in the decay width is: 


( Ar )cW. = -i^~9o,oo( r = 0)D 0 2 0i+ _b + _, + _^_(0) 

y f, ^+-,n(0) 

X m n (E 0 - E n ) ' 


(5.1) 


Using the nonrelativistic formulas 


E n = 


m n a 


mta 


4(n+l) 2 ’ 8jr („' + l) 3 


(5.2) 


one finds: 

q, _ q, 

( A -^)djscr. = — 4^ m7r ^+--+-^° = g - (^ a 0 + a o)r(b (5.3) 

where in the last expression we have neglected the pion mass difference and introduced 
the strong interaction scattering lengths (Q. 


5.2 Contribution of the continuous spectrum 

For the evaluation of the contribution of the continuum states, we introduce the nonrel¬ 
ativistic wavevector modulus k = \/Em n and normalize wave functions as: 

J d 3 x^,(x)'0 fc (x) = 2nS(k — k'). (5.4) 


The shift in the decay width is then: 

( Ar U, = -2r 0 F + _ i+ _ 


with 146 


hM0)l 2 = 


r dk |^fc(0)| 2 
J 2n {k 2 + m^a; 2 /4) ’ 

A-Kin^ak 


(5.5) 


(5.6) 


47 t (1 — exp(—7 rm n a/k)) 

It is easily seen that ( AT) cont diverges linearly in the ultra-violet region. In dimen¬ 
sional regularization, linear divergences being equivalent to zero, the divergence that 
survives in Eq. (|5.5|) is logarithmic: 

k . 


(AT) 


cont. 


on T7 Tn n cy . 

— 2r 0 U— ln(- 


(5.7) 
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The origin of this divergence is related to the singularity of the three-dimensional delta- 
function that characterizes the potentials in x-space. In quantum mechanics, potentials 
that are proportional to three-dimensional delta-functions must be regularized through a 


renormalization of the coupling constant 47| or a self-adjoint extension of the hamiltonian 


48|. In both cases, an additional unknown parameter appears in the spectrum. On the 


other hand, a bound state equation, based on field theory, should not introduce new 
unknown parameters into the set of parameters already defined or fixed by the field 
theory itself. This is why one should expect the cancellation of the divergence by 
some other terms present in the potential. 

The constraint diagrams corresponding to the process 7r + 7r“ —» 7r°7r° with one photon 
exchange and two loops (Fig. 2) play this role. They have overall logarithmic divergences 
that cancel the one appearing in (AT) conf . Therefore, these two types of contribution 
should be considered together. 

We now consider the contributions of the above constraint diagrams. At the vertex 
— 7T°7T°, the amplitude that contributes is A(s\t,u) [Eqs. (|3.5| ) and (|1.1| )1; in lowest 
order, it depends only on the variable s and therefore can be factored out of the integrals. 
At the vertex 7r + 7r“ — 7r + 7r“, it is the amplitude A(s\t,u)+A(t,\s,u) that contributes. 
However, the contribution of the variable t, present in A(t\s,u), is of order 0(a 2 ) when 
considered in the bound state domain (cf. also Appendix A); therefore, it can be omitted 
and the amplitude A(s\t, u)+A(t\s, u) can also be factored out. We use for the potentials 
arising from the constraint diagrams the notation C referring to “constraint”, 

n to the number of loops, m to the number of exchanged photons and p to the number of 
the constraint factors go. 

We first consider the constraint diagram of Fig. 2-c. The corresponding potential is 
(in the Feynman gauge): 

d 3 k T cl 3 k' T 




C( 2,1,2) 


00 ,+- 


= e 


2y/s 


Voa,+-V 


x 


(27t) 3 (27t) 3 
(2 pi — k T — k' T ).(2p 2 + k T + k’ T ) 


—i 


{{pi — k ,T ) 2 — ml + ie)((pi — k T ) 2 — ml + ie) ( k T — k ' T ) 2 + ie 
d 3 k T 1 


-ie 2 Vnn. 4 —. 4 — 


(27t) 3 (pi — k T ) 2 — ml + ie 
x|(4pi.p 2 + 4p T .k T )F (1)c (p ll k T ) + 4 (p T - k T YF [3)C (p 1 , k T ) - A c (pi)|, 

(5.8) 


15 















with: 


FWC ^ = 2hi 


d 3 k' T 


2^/s J (2n) 3 ((pi — k ' T ) 2 — ml + ie)((k T — k ,T ) 2 + ie) 
F (1)c (pi -fc T ,0) 
i 1 2 


arctan, 


4^471 ^-{ p t - AT) 2 \ - & o(s) ’ 


— (p T — k T ) 2 


(5.9) 


- 6=J 


d 3 k' T 


h’ T 

V 


2^/s J (2n) 3 ((pi — fc ,r ) 2 — m 2 + ie)((k T — k ,T ) 2 + ie) 
klF ( » c ( Pl - k T , 0) + Fj 1)c ( Pl - fc T .0), 


Jf )C (Pi,0) = 


=rC 


' Z ' rp I S 

8~s Pfl \lv 


3 uT 


62 arctan W 


b 2 


bl 


6 2 


#"(Pi) = 




- 68 . 


(5.10) 


(5.11) 


2yT J (27t) 3 (pi — k T ) 2 — m 2 + ie 2yT 4n 
[The integrals are finite in dimensional regularization. We have kept in FjA c the leading 
terms only.] 

In Eq. ( |5.8| ), it is only the term proportional to pi.p 2 , through PilP 2 L, that gives the 
surviving 0(a) contribution. This observation makes also clear that the consideration of 
other covariant gauges would also yield the same leading term. Therefore, the subsequent 
results are gauge independent. One obtains: 

d 3 k T 1 1 


Kn ( l’±’ 2) = -m n aV { 3 o,+-E + _ >+ _ 


arctan 


-AT 2 


(5.12) 


“, + _ V -6 1 

We next consider the two constraint diagrams a and b of Fig. 2. Since the amplitude 
A(s\t,u)+A(t\s,u) has been factored out of the integrals, the two diagrams give equal 
contributions. The four-dimensional loop corresponds to the electromagnetic radiative 
correction of the scalar vertex. Its expression will be calculated in Sec. 6 (and Appendix 
D). It contains, among other terms, a term that is the opposite of the function F^ c 
[Eqs. (|5.8|) - (|5.9| )1 we calculated above. [It is the first term in the expression of F T in 
Eq. ( p . 1 1| ) , where b 2 = —p T2 should now be replaced by — (p T — AT) 2 .] Its other terms 
contribute to nonleading effects in a. 

Therefore, the contribution resulting from the sum of the two diagrams (a) and (b) of 
Fig. 2 is: 

V$f M) = -2v£jy>. (5.13) 


16 
































The contribution to the shift in the decay width resulting from the three constraint 
diagrams of Fig. 2 is then: 

d 3 k 1 


(AT) c(2 ’ lll+2) = 2m n aT 0 V + . )+ . 


1 , 2k 

arctan (- 


m^a 


(5.14) 


(27t) 3 (k 2 + m\a 2 j 4) k 

where we have replaced —5 q = ml — p\ L by its eigenvalue m^a 2 / 4. 

(Ar) c has an ultra-violet divergence that is the opposite of that of (AT) cont [Eq. 
(fTTPI. Therefore, the sum (AT) C '+(Ar) corit is hnite. One obtains: 

(AT) c(2 ’ 1>1+2) + (AT) — 1 


' cont. 


2m n aT 0 V — 


p OO 

x / dx 
Jo 


X 


X 2 + 1 


2 x 

— arctan a; H- 

7T 7T 1 


p—2ix/x 


a „ — 


-1.7m n — T 0 E + + _ = l.la(2ao + ao)T 0 . 

Z7T 


(5.15) 


[We have subtracted from (AT) corit its linear divergence, which is null in dimensional 
regularization.] The constraint diagrams have played the role of an effective cut-off of the 


divergence of (AT) cont . Had we cut the integral in Eq. (|5.5|) at the value k = m n a , we 
would have found a result close to that of Eq. (|5.15|) . 

In the previous cancellation mechanism of divergences we used dimensional regular¬ 
ization and hence ignored linear divergences. The latter are present in conventional cal¬ 
culations of integrals and for such type of calculation they should be taken into account. 
The point is that in this case the constraint diagrams themselves are linearly divergent 
and these divergences should be isolated and grouped essentially with those coming from 
second-order perturbation theory. It is not difficult to show that the linear divergences 
cancel out among themselves and the physical results are not affected by their presence. 
Details of the derivation are presented in Appendix C. The cancellation mechanism of 
divergences is therefore regularization scheme independent. 

The total amount of the strong interaction corrections is given by the sum of (AT) rfiscr 
[Eq. (p|)] and (AT) C(2 ’ 1,1+2) + (AT) conf [Eq. (|5T5|) ]. One obtains: 


(Ar ), tr . = (Ar) iiiCT + ((Ar) co „, + (Ar) c(2J - 1+2 >) 

= 1.5a(2ag + ao)T 0 . 


(5.16) 


The value of (AT) str depends on the values of the scattering lengths and hence also of 
the quark condensate parameter 2mB 0 [Eq. (|1.1| )]. In the standard case ( 2rhB 0 ~m^) 0, 
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g one has 1.5(2ao+Oo)—0.55. In the other extreme case (2rhB 0 = 0), one has 1.5(2ao+do) 
Therefore, in all cases the correction (AT/T 0 ) str .. is bounded by the values: 


~0.9 22 



0.55a = 0.004 
0.9a = 0.0065 


(2 rhB 0 = ml), 
(2mB 0 = 0). 


(5.17) 


6 Electromagnetic radiative corrections 


In this section we calculate the electromagnetic radiative corrections arising fom pion- 
photon interaction and contributing to the amplitude A4 em ' of the decomposition ( |4.3| ). 

We hrst evaluate the unrenormalized quantities. The lagrangian corresponding to 
the pion-photon interaction is obtained by appropriately incorporating in the external 
vector current of the chiral effective lagrangian the photon held. Using for the held U the 
representation (|3.4| ), the corresponding lagrangian becomes in lowest order: 

C = d^n 0 + (dp + ieA fl )n + (8 fl - ieA^)n~ 

~\ m l ( 7r ° 7r ° + 2vr + 7r-) - {8.A) 2 

+ ^‘( 7r ° 2 + 27r + 7T")cI / "(7r 02 + 2 ' k + h ~) 

- —2 (2mi?o)(vr 04 + 47r + 7r _ 7r°7r° + 47r +2 7r -2 ), (6.1) 

8 F 0 

2mBo being dehned in Eq. (|1.1|) . 

We begin with the self-energy correction, which is (photon tadpole contributions are 
null in dimensional regularization): 


- «E(p) 


r d 4 k ((2p — k) 2 — £\(2p — k).k] 2 /k 2 ) 

I (2i r) 4 ((p — k ) 2 — ml + ie) ( k 2 + ie) 

4 p 2 R {1 \p) + (1 - £)T(P) - 4(1 - 0P^\P) 


4 &YR$(p) , 

( 6 . 2 ) 


where ^ represents the gauge parameter. The definitions and expressions of the functions 
R and T are presented in Appendix D. The expression of the unrenormalized self-energy 
is given in Eq. ( p.8| ). The unrenormalized Green’s function of the charged pseudoscalar 
density becomes: 


G+(p) 


pZ 


t(2B 0 F 0 ) 2 

(ml + (AmJ) n ) + ie 


(1-A Z + (p)), 


(6.3) 
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with: 


( Am JW = (^ + ^K> 

— AZ + (p) = (2 + ^(l + (1 + |) ln(^) + '(I + |)^l 


(6.4) 

(6.5) 


2'~~' K \ 2J 2 2 V A 2 

£ and A 2 being defined in Eqs. ( p.6|) and (|D.7|) , respectively. 

We next pass to the vertex function (Fig. 3a). Omitting from its definition the 
7 r + 7r - —> 7r°7r° amplitude A(s\t,u), which at the tree level does not depend on t and u 

[Eq- (0)]> it is: 

d A k ((2pi - k).{2p 2 + k) - £(2pi - k).k{2p 2 + k).k/k 2 ) 


A(pi,p 2 ) = ie 


(27t) 4 ((p i — k) 2 — m 2 + ie)((p 2 + k) 2 — m 2 + ie)(k 2 + ie) 
ie 2 [-R (1 \pi) + (Api.p 2 + p 2 2 - m 2 )F (1) (pi,p 2 ,0) + 2p) t iA 1) (pi,p 2 , 0) 
c) 

- Kv-v I n, ) -I- 1 777 * — nt KZnV ^ 


-*e 2 £ 




i? (1) (pi) + (m 2 — pl)2piFff\pi,p 2 , 0) 


dpi 

-2p^F^\ Pll p 2 ,0) -R( Pl ,p 2 ) , (6.6) 

with the definitions and expressions of the functions R , A and F given in Appendix D. 

To this order of approximation, the strong vertex part of A is a scalar vertex and 
therefore A satisfies an identity typical of mass operators, obtained by differentiation 
with respect to m 2 P5fl . It reads: 

<9E( Pl ) 


HPu -Pi) = 


dm 2 


(6.7) 


The unrenormalized A at leading orders is given in Eq. ra¬ 
in passing to the scattering amplitude, one must associate with the vertex function 
the wave function renormalization factors of the external particles (the square-root of 
(1 — A Z + ) of Eq. (|6.3| )). This leads to the combination: 

aAp 1 .p 2 


A(pi,k) - Iaz+m - 1a Z+(p 2 ) = (F _ !_) + 2 


— arctan - 


I b 2 

-N' 


( 6 , 8 ) 


[We have kept dominant terms only.] Notice that the gauge parameter and the logarithmic 
functions have been cancelled out. 

With the above quantity we must associate the contribution of the constraint diagram 
arising from the interference term between the strong interaction amplitude and one- 
photon exchange diagram (Fig. 3b): 

i f d 3 k T ((2pi — k T ).{2p 2 + k T ) — £(2pi — k T ).k T {2p 2 + k T )/k T2 ) 


R c {p\pp 2 ) = ie' 


2yfiJ (2vr) 3 


((pi — k T ) 2 — m 2 + ie)(k T2 + ie) 


19 



















-ie 2 i 


R°(pi ) + 4pi. P 2F wc ( Pl , 0) + Ap^F^iPu 0)' 
2X 2 p T,x F^ c (p 1 , 0) + 20) - Tf ( Pl ) 


(6.9) 


the definitions and expressions of the functions R and F c being given in Eqs. (|5.9|) - (|5.11| ) 
and Appendix D. One finds: 


A c (pi,p 2 ) = “f ^^/Jarctan ] j-^. (6.10) 

Adding this term to the quantity in Eq. (|6.8|) yields the net effect of the unrenormalized 
radiative corrections on the potential Voo,+- (relative to 1): 

1 1 py 

A(P1,K) - ? AZ + (Pi) - -AZ+fe) + A c ( PllP2 ) = (-—), (6,11) 

2 2 el 7T 

which is a gauge invariant quantity. On the other hand, the constraint diagram has 
cancelled the term proportional to the arctan function present in A, which otherwise 
would give a contribution of order O(a 0 ). 

We next consider the contributions of the 0(e 2 p 2 ) counterterms present in the chiral 
effective lagrangian. The complete expression of the corresponding SU(3) x SU( 3) la- 
grangian in the Feynman gauge (in the standard scheme, to which we stick henceforth) 
is given in Ref. [£h|, where the coefficients of the various terms are designated by Ki 
(i — 1,..., 17). These terms are of two categories. Those which are counterterms to the 
0(e 2 p°) lagrangian ( |4.1| ) and the infinite parts of which are proportional to C, and those 
which are counterterms to the radiative corrections of the pion-photon interaction. Since 
the (0(e 2 p°) term is included in the lagrangian contributing to the amplitude J\A str ' +97 , 
so must be done with the corresponding counterterms. Therefore we have to select from 
the 0(e 2 p 2 ) lagrangian those terms the infinite parts of which are not proportional to the 
constant C. 

The 0(e 2 p 2 ) lagrangian in the SU{2) x SU{2) case was recently presented in Ref. 
||5C)| (in the Feynman gauge). Since we are treating in the present work the SU(2) case, 
we use for our subsequent calculations the latter lagrangian, indicating, when necessary, 
the correspondence between the SU(2 ) and SU{3) coefficients. The relevant part of this 
lagrangian for our problem is (in standard notations): 


£(eV) = F 2 jfc 2 < QUQU ] >< dYdY > 

+fc 3 ( < U ] d„UQ >< U'dPUQ > + < >< d^UU^Q > ) 
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+k 4 < U ] d lt UQ >< d p UU ] Q > +k 7 < QUQU ] >< X U ] + U\ ] > 

+k 8 < (U ] x. ~ x ] U){U ] QUQ - QU^QU) > 

+k 10 < Q 2 >< duU&tf > +k u < Q 2 >< XU ] + Ux 1 > ], (6.12) 

where Q is the quark charge matrix. The coefficients k have the following decompositions: 


where A has a pole in d — 4 dimensions, 


A = 


1 


A + fcJXh), 

(6.13) 

^(ln(4vr) + T'(1) + l)) , 

(6.14) 

fiq = — 2Z, k 7 = j + 2Z, 

K - 1 z 7 - c 

K 11 - “4 _ 5’ Z - F*- 

(6.15) 


167T 2 \d — 4 

and the coefficients k have the values^: 
k 2 = 2 Z, k 3 = -|, 

— g — Z, Kio ~ ~ 20 — ¥’ 

The /^-dependence of the hnite part k [ (/t) of k t is fixed by the prescription that its sum with 
the ln(/i -2 ) term of the corresponding loop diagrams be /^-independent; the corresponding 
multiplicative coefficient is /q/(327T 2 ). The lagrangian ( |6.12|) , together with the coefficients 
( |6.150 , is associated with the SU( 2) x SU(2) version of the chiral effective lagrangian of 
Ref. in which the multiplicative factors have been arranged so that the coefficients U 
keep the same values as those of the 0(4) version [Eqs. (|7T Q -(|7.3 Q 1. 

The lagrangian (|6 .1 2Q takes the following form in the representation Q3.4|) , keeping 
only pion helds that contribute to the scattering problem to the present approximation: 

10, fc. 


£ {ep) - e 2 [k 2 + y(y§ + M) 


d u 7T 0 d p 7i° + 2d u n + d fJ 'n~ 


+jr 2 {\ay 2 + ((p'v 02 + d^+n-)) 


4e 


T k 2 TT + 'K (d^lT°d^TT° + 2<9 m 7T + <9 / V ) 


^0 

2 2 
+ -^(2 k 3 - k 4 )(n + dp n~) - -^(2 k 3 + fc 4 ) 


o 


F 2 d^°d^° 


1 


—n [>2 d l _ l 7r^d tJ '7T [> — 2 tt + 7t d^d^i r u + ^<9 /t 7r 02 <9 M 7r 02 + <9 At (7r 02 )<9 M (7r + 7r ) 


5e 2 

- —(fc 7 + k u )8rhB 0 


-7T° 2 + 7T + 7T H-^(7T° 2 + 27T + 7T ) 2 

2 8Fq v 


—~p2 (^(^ 7 — 2fc s) + ^(^7 + 2fc 8 ) )8mB 0 


Fq TT^TT — 7T + 7T (-7T° 2 + 7r+7r ) 


(6.16) 


4 We thank M. Knecht for pointing out to us the presence of the factor —Z in k&. 
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The coefficients (k 2 / 10+fcio), k 8 and ( k 7 +2k 8 ) contribute, at least concerning their infinite 
parts, only to the pion-photon interaction effects, while the coefficients k 2 , k 4 , ( k 7 + ku) 
and (kj — 2kg) contribute to the quark-photon interaction effects (insertions of the mass 
term (|4.1|) ). We shall admit that this separation also holds for their finite parts. 

The 0(e 2 p 2 ) lagrangian (|6.12| ) introduces new terms into the relationships between 
the pseudoscalar densities and the pion fields. These are: 

P a = 2B 0 F 0 (l + — (k 7 + k u ) + 4e 2 k 8 5 a± -2e 2 ((k 7 -2k 8 ) + 3(k 7 + 2k 8 ))—^)n a , (6.17) 

where a = 0, +, —. The three-pion term yields additional diagrams for the scattering 
amplitude, which, however, contribute only off the mass shell. The term proportional to 
the coefficient k 8 alone contributes to both the pion-photon and quark-photon interaction 
effects. For this reason we shall separate this coefficient into two parts by writing, with 
an obvious notation, k 8 = fcg 7 + fcf 7 . Although this separation is ambiguous for the 
corresponding finite parts, it will not show up in physical quantities. 

To obtain the renormalized Green’s functions of the pseudoscalar densities due to 
pion-photon interaction, we select among the above expreesions the quantities propor¬ 
tional to the relevant coefficients. Furthermore, since we are considering here on-mass 


shell expressions, the three-pion term of Eq. (|6.17|) can be discarded. One finds (in the 
Feynman gauge) the expressions: 

i(2B 0 F 0 ) 2 (l - A Z°) 


G U (P ) = 


(Am: 


-A Z° 

2 \ 


p 2 - (ml + (Am( 0 )i 7 ) + ie' 

+ iofcy + ie%, 


/7T7 


= 


g - (^2 + lO^io) + 4e 2 &3 


(6.18) 


G + (p ) 

-AZ + 


(Am) 


■+ ) 7I"7 


i(2B 0 F 0 ) 2 (l- AZ+) 


ml + (Amjt) n ) Fie 
2 ^ \2 ^2 


a « , /ml, a .ml a / ,mt. \ 

+ lOfeJo) + 8e 2 C', 


= 


S “ £( + + 10l “ ) + ^ + 2k *\ 


(6.19) 
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The difference between the charged and neutral pion masses due to pion-photon in¬ 
teraction is: 


(A ml) n -f = (AmJ+Jry — (Arajo) 


7T7 


7a 3 a/, ,m 2 


4 „(M^) + l)-4An6e 2 ffi + 2*:) 

To obtain the renormalized scattering amplitude of the process 7r + 7r - - 




( 6 . 20 ) 


7r°7r°, we con¬ 


sider the sum of the unrenormalized four-point vertex function and of the corresponding 
constraint diagram contribution: 


W>+- + 


= (2B 0 F 0 )- A (l+Ae 2 kr)~ 2 A(s\t,u)(l + A + A c ), (6.21) 


J 7T7 


where A(s\t,u) is defined in Eq. (0), A in Eq. ( p,19| ) (taken here in the Feynman 
gauge) and A c in Eq. (|6.10| ); the other multiplicative factors come from the differences 
between the pseudoscalar densities and the pion fields [Eqs. (|6.17|) ] . The lagrangian ( |6.16|) 
provides the following counterterms: 

.2 rt2 


(W / oo,+-)tt 7 — (25oFo) (-(^2 + lOfcio) — 8fc 3 )s 

r n ' 


'9 


3 

+8(fc 3 + ~(k 7 + 2kg)^)2mB 0 


( 6 . 22 ) 


These remove the divergences of the unrenormalized vertex (|6.21|) . The sum of the two 
quantities is then multiplied by the renormalized wave function renormalization factors 
(2£? 0 Eo) 4 (l — AZ°/2) 2 (1 — AZ+/2) 2 [Eqs. (|6.18|) - (|6.19|) 1. We also incorporate in the 
mass term of the lowest-order amplitude (|1 . 1|) the neutral pion mass shift (|6.18|) . The 
renormalized scattering amplitude, together with the constraint diagram contribution, is 
then: 


Ad 00,4— + ■Moo'+d ^ 


(s - m 2 0 ) 


J 7T7 




a 3a / ,m 2 \ ■>, 

1 - 7 — - — ln(—+ 1 - 4 e 2 k 
2t r 4 tt V v /i 2 ’ > 


2 A(kl + Wk\ 0 ) 


o 

r,mz 

+ 2e 2 -| 

^0 L 


2e 2 k r 3 + 3e 2 (k r 7 + 2k r 8 ) 


(6.23) 


The corresponding modification of the decay width is: 


(A T)^ = 


- 7 - ^ (In(^) +1) - ^( k r 2 + 10A* 0 ) - ^fk r 3 + 4 e 2 (k r 7 + 2kl) 


T 0 . (6.24) 
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To estimate the correction (|6.24| ) one needs to know the values of the coefficients 
k r . At this time they are not yet well established, but several evaluations with different 
methods exist and can be used to have an idea of the order of magnitude and sign of the 
correction. The correspondence of the above coefficients k with the SU(3) coefficients K 
is (valid also for the finite parts): 


+ 10 /cio — 12 A i + IOA5, k 3 — A3, 

-(/C7 + 2k$) = — -A s + ^10 + -^n- 


(6.25) 


In ref |51|, the K r, s have been evaluated using a resonance model for the saturation 
of sum rules. The corresponding values are in units of 10~ 3 and at the scale p = m p : 
K[ = -6.4, iq = 6.4, Kl = 19.9, iq = iq = AT[ 0 = 0, AT n = 0.6. 

In Ref. [[5^] values of some of the A' r, s have been evaluated with the aid of the 1/N C - 
expansion method, the Extended Nambu Jena- Lasinio model and perturbative QCD. 


The results are not completely compatible with those of Ref. |5T[ . In particular it is found 
(at similar scales and in units of lO” 3 ): Kg = —1.3, K\ a = 4.0; the other results involve 


K r, s with even indices. Estimates of K\ 0 and K[ x have been presented in Ref. |53| on the 
basis of sum rule analysis. The value of K\ 0 found there is compatible with that of Ref. 
f>2| , while for K\ x the value 2.9 x 10 -3 is found (with the choices p = po = mv = wa/ \/2). 


The latter values, which affect the combination (&£ + 2fcg), do not seem, however, to 
qualitatively change the estimate of the correction (|6.24|) evaluated with the values of 
the coefficients of Ref. 0- We shall therefore present the numerical estimates with the 


values supplied by Ref. [51 


We find (with the mass scale p = m p ): 

(Am f o) x1 ~ —0.01 MeV, (Am 7r +) 7r7 ~ 0.43 MeV, (A m n ) ni = 0.44 MeV ., 

(AT) ¥7 = —0.20aT 0 = —0.0015T 0 . (6.26) 

We observe that the pion mass shift is of the order of 10% of the observed mass shift 
A= 4.6 MeV and has the correct order of magnitude for an 0(e 2 p 2 ) effect. The 
decay width correction is found almost equal to zero. The counterterm coefficients k r of 
Eq. (|6.24|) give a negative contribution and thus have the tendancy to lower the value 


of (AT) W7 ; the contribution of the Erst two terms of the right-hand side of Eq. ( |6.24|) is 
0.0025T 0 . 
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7 Electromagnetic mass shift corrections 

This section is devoted to the evaluation of the electromagnetic mass shift corrections 
contained in the amplitude 7leA4 s 0 t Q^T / [Eqs. (f4~7|)- (|4.8| )1. We shall proceed in two steps. 
First, we shall determine the effect of the insertions of the mass term (|4.1|) in internal 
charged pion propagators. This will bring us to the strong interaction amplitude cal¬ 
culated with individual momenta fixed at the charged pion mass, but with the mass 
parameter fixed at the neutral pion mass. In the second step, we shall determine the dif¬ 
ference between the latter amplitude and the one calculated with individual momenta and 
masses fixed at the charged pion mass, with which the numerical values of the scattering 


lengths are calculated in the literature ||15| , |23 


For the evaluation of the effect of the insertions of the mass shift term ( |4.1| ) the coun¬ 
terterm lagrangians needed here are the quark-photon part of the lagrangian (p . 1 2 | ) and 
the standard 0(p 4 ) lagrangian part of the latter survives in the mass shift countert¬ 
erms. We therefore begin by writing the relevant part of this lagrangian for our problem 
(with standard notations): 


k 


U 


= j < d^d p U > 2 +j- < dprfdJJ >< d p U ] d v U > 

+ ii < x ' u+xU] >2 + I < d ^ d " U + dpXd " U] > 

The coefficients l have the following decompositions: 


(7.1) 


h — 7 A + 


(7.2) 


where A is defined in Eq. (|6.14|) and the coefficients 7 have the values: 


12 1 

71 = o’ 72 = g, 73 = 74 = 2. 

For the scattering problem the above lagrangian becomes: 

C^ = ■^(d fl 'K°d p n 0 + 2d fl 'K + d p 'TT~)(d u 7[°d u 'iT 0 + 2d u TT + d u n~) 


(7.3) 


+ -^(d li TT°d u 'K° + d f ,TT + d u 'K +^7T d v TT + ) 


X^TT^TT 0 + d^TT+ff'TT- + d P Tl~d U Tl + ) 


u 


-^(2mi?o) 2 (vr 02 + 27r + 7r-). 
-^0 


(7.4) 
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This lagrangian provides the following relationship between the pseudoscalar densities 
and the pion fields: 


/ 2 h 

P a = (2B 0 F 0 ){l + ^2mB 0 (l 


vr- , 

2Fl' 




(7.5) 


The three-pion term does not contribute to the present appreximation to the two-point 
Green’s function renormalization, but provides an additional contribution (a three-propa¬ 
gator term) to the four-point Green’s function. Also, the term proportional to Z 3 in 
Eq. (|7.1|) provides in addition a two-propagator term to the four-point Green’s function. 
These contributions should be taken into account for the off-mass shell expression of the 
scattering amplitude. 


The finite parts of wave function and mass renormalizations proportional to C/Fq 


111 


notations similar to those of Eqs. ( f3.18| )-( |0^ )) are: 
—2e 2 


A Z° 
.2 




k r 2 


-2 e 2 m 2 


K 


20 

— (A^ + ArJi) 


fw +k 


ii; 


(7.6) 


— AZ + = ~2e 2 




4 k™ 


2 C„ 

774 4 
^0 - 


(Am 7r +) q7 — 


2e 2£L _ 


2 C 1 

Ff 16^ 


777 ^ 1 0 

ln(^) + k r 2 --(k r 7 + k[ 1 )-(k 7 -2k r 8 ) 


(7.7) 


(We have also incorporated in (Am^+) g7 the lowest order term ([0])). The difference 
between the charged and neutral pion masses due to quark-photon interaction is: 


(A ml) qi = (A ml+) qi - (Am 2 0 ) 


97 


9 C 

= 2e V 

-^0 


2e 2 m 2 


C 1 

Fg 167 t 4 v ' ji 


(31n(^) + l )+k r 4 -(k r 7 


2 kl 


(7.8) 


To evaluate the change in the scattering amplitude due to the pion mass shift (|4.2| ), 
we first calculate the pion loop contributions with the physical pion masses in the internal 
propagators. There are diagrams with tadpole loops at the four-pion vertex and four di¬ 
agrams with loops with two-pion propagators, one with two charged pions, one with two 
neutral pions and two with one charged and one neutral pions. To the contributions of 
these diagrams one adds the remaining contributions of the pseudoscalar densities to the 
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four-point vertex function [Eqs. ( |S.17| ) and ([775])]. The divergences of the above contribu¬ 
tions are cancelled by those of the counterterm lagrangians ( |6.16| ) (the part proportional 
to the factor C/F$) and ( |7.4|) . One multiplies the result by the wave function renormal¬ 
ization factors (1 — AZ°/2 ) 2 and (1 — AZ + /2) 2 coming from the contributions ( |7.6|) - (|7.7|) 
and from those of the strong interaction limit § and isolates in the resulting expression 
the part proportional to C/Fq. One also incorporates in the mass term of the lowest-order 
amplitude (| 1 . 1 |) the neutral pion mass shift ( |7.6| ) . One finds in the limits t — u — 0 (with 
the notations m^ = m^o and Am 2 = (Ara 2 ) g7 , Eq. (|7.8| )): 


F 


t (A77eAf 0 o, + -) 


97 


( s - ml)(s + pI+P 2 2 - 2 ml) + 77e(Q + (a) - Q 00 (s 

+Y^( s - m l)( 2 - ne Q 00 ( s )) Am l~-^(pi+p 2 2 )(pl+pl)^^- 
+2e 2 F 2 (3k r 2 - 2 k\)s - 2e 2 F 2 {2k r 2 - k r 4 )(pl + p\) 

*■2e 2 F 2 (k r 2 - kl - (k r 7 - 2kl))ml - 2e 2 F 2 (2k r 2 - k r 4 )(s - m 2 ) 

+2 e 2 F 2 (k r 7 - 2k r s )(p 2 3 + p\- 2m 2 ) 

- 4 / 3(73 + Pi - ml)Ami ~ 2l l( s - ml) Ami. 


32tt 2 


(7.9) 


(Details of the calculations can be found in Appendix E. Q(s) is defined in Eq. ( |E.3| ). 
Only relevant hrst-order terms in Am 2 have been kept.) 

With the mass shell conditions p 2 = m 2 + (i = 1,... ,4) and s = 4m 2 + , the above 
expression becomes, to hrst-order in Am 2 : 


F 0 4 (A77eM 0 o,+-(s = 4m 2 + )) 


97 


11m 2 

167T 2 


■A/n 2 - ( 4/3 + 6 / 4 )m 2 Am 2 


+2e 2 F 2 ((^ - 2 k r A ) + (k r 7 - 2kl))ml. (7.10) 


The relationships between the k\'s and the 5/7(3) coefficients Kl [|45|] are (m# is the 
kaon mass): 


kl = Kl 


•2 - “2 + ^-gTj(ln(^) + l), K = -K, 


1 Q / 777 ^ \ 

*? + * [1 = P6K; + 5K; + 5Kl 0 )-— 2 (ln(J-) + l), 


kl — 2kl = 10Kl - 
7 8 8 16tt 2 


3 /, ,m 2 


C“(f) + i)- 


(7.11) 
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Numerical values for the coefficients K% are presented in Ref. |)]J. They are in units 
of ICR 3 and at the scale fi = m p : K £ = —3.1, K r A = —6.2, A'g = 8.6, Kg = Kg = K r m = 0. 
The values of the coefficients l\ can be found in Ref. (15J. They are in units of 10“ 3 and 
at the scale \x = m p \ l\ = —5.4, l 2 = 5.67, l 2j = 0.82, l r A = 5.6. One finds the following 
corrections in the scattering lengths and the decay width: 

(Ar)„ 


(A(oq - «/ ( ) ),, = 0.0005, 


= 0.0035. 


(7.12) 


Once the correction due to the mass shift ( |4.2|) is isolated, the amplitude IZeKigg + 'A 1 
reduces to the strong interaction amplitude calculated with individual momenta fixed 
at the charged pion mass [Eq. (|4.8|) 1 and the mass parameter fixed at the neutral pion 
mass. Numerical values of the strong interaction scattering lengths are calculated in the 
literature with the individual momenta and the mass parameter fixed at the charged 
pion mass [|T5| [23|| . It is therefore necessary to evaluate the difference between these two 
amplitudes. It can be obtained from the 0(p 4 ) off-mass shell expression of the scattering 
amplitude given in Ref. H (we neglect for the moment the mass shift effect in the 0(p 6 ) 
term). Replacing in the 0(p 2 ) term (|1 . 1|) the mass parameter in terms of m 2 0 with the 
use of the relation [01: 




(7.13) 


one finds for the shift at threshold from the amplitude calculated with the charged pion 
mass: 


(ATleMoo, + -(s = 4 m£+)) 


Am-n 


Ami 


1 + 


mz 




(- 


9 


32vr 2 
+4/3 + 12 ll ) 


11 ,771. 


(7.14) 


Designating by _the strong interaction amplitude (without mass shifts), one 


has at threshold [Eqs. (|3T5|) , ( fTfip and ([37^)1 : 


327T 

A45J+-( S = im l) = - 5 -K - 


(7.15) 


where the a’s are the scattering lengths calculated up to two loops |15|, Q of the chiral 
effective lagrangian. In the standard scheme one has the value [[|, [L5|] : (a° — ao) = 0.258, 
obtained with charged pion masses. Eq. ([7.14 ) yields the following corrections for the 
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scattering lengths and the decay width (we take F 0 = 88 MeV, m w + = 139.57 MeV, 
Amj = 4.6 MeV): 


(A(ag - ao)) Am ^ = 0.0083, 


(AT) 


Am-n 


= 0.064. 


(7.16) 


The contribution of the 0(p 4 ) term represents 40% of this correction, indicating an in¬ 
crease of its relative strength by a factor of 2 with respect to the corresponding situation 
in the amplitude; this might be understood as a consequence of the increase of the powers 
of mass and momemtum terms in higher-order terms. In this case one should expect a 
correction of the order of 3 x 5% coming from the 0(p 6 ) term, which would bring the 
correction (|7.16|) in the scattering lengths to 0.0095. 

The effective scattering length corresponding to the amplitude IZeMw'+V ( |4.7| ) is 
then (without the use of the 0(e 2 p 4 ) correction coming from the 0(p 6 ) term): 


K a o)str.+q'Y — (®o a 0 )str.(n+) T (A(Oq a o))qi (A(®o a 0 ))Ami — 0.267. (7.17) 

(The inclusion of the estimated 0(e 2 p 4 ) correction coming from the 0(p 6 ) term would 
bring this number to 0.268.) 

It is also possible to extract the value of the constant C from the pion mass shift 
formula (|7.8|) . Taking into account the mass shifts due to pion-photon interaction (0.4 
MeV, Eqs. ( |6.20| ) and ( |6.26|) ) and to isospin breaking (0.2 MeV [[|, |J) and the observed 
pion mass difference (4.6 MeV), one finds (Am^) qi ~ 4.0 MeV, leading to: 


C = 4.1 x 10' 5 GeV 4 . 


(7.18) 


This value is in agreement with the central value found in Ref. [Q at a slightly higher 
mass scale ((4.2 ± 1.5) x 10^ 5 GeV 4 , at p = 0.85 GeV). The corresponding O(e 2 p 0 ) mass 
shift is (A m n )^ p °' > ~ 3.6 MeV, the remaining 0.4 MeV being produced by the 0(e 2 p 2 ) 
correction [Eq. ([7.81)1 , with (Am jr o) g7 = —0.13 MeV. 


8 Summary 


The total amount of 0(a) corrections to the lowest-order formula of the pionium decay 
width can be represented in the following form: 


T 



Am,)(i + (AEA)(i + AIA)(i + AEkm + (Ar)aro - ), 

2m,r+ v r 0 yv r 0 > 


( 8 . 1 ) 
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where T 0 is the lowest-order decay width, Eq. with the strong interaction scattering 

lengths calculated up to two-loop order of the chiral effective lagrangian with charged 
pion masses; (AT) s j r . is the correction arising from the second-order perturbation theory 
contribution of the interaction potential; (Ar) T7 arises from the radiative corrections due 
to pion-photon interaction; (Ar) g7 is the correction coming from the pion mass shift due to 
quark-photon interaction; (Ar)/\ mir is the correction coming from the shift in the strong 
interaction amplitude due to the passage from the charged pion masses to the neutral 
pion masses with individual momenta fixed at the charged pion mass. The corresponding 
numerical values are (in the standard scheme): 

~ 0.004, _ _0.001, 

To 

~ 0.003, ^ r ) Am ' ~ 0.064. (8.2) 

To 

These numbers are subject to uncertainties coming from the numerical values of the 
coefficints and l\ of the counterterm lagrangians and from the higher loop effects not 
considered in the pion-photon and quark-photon interactions, corresponding to 0(e 2 p 4 ) 
terms. These corrections should not qualitatively change the above numerical estimates. 
One finds for the lifetime the value: 


(AT) s t r . 

r 0 

(Ar)„ 

U 


r = 2.97 x 1(T 15 s, 


(8.3) 


to be compared with the value r 0 
a§) = 0.258. 


3.19 x 10“ 15 


found from formula (|1.2| ) using 


o 

o 


The evaluation of the corrections in the generalized scheme is more complicated be¬ 
cause of the presence of additional terms in the counterterm lagrangians. However, their 
order of magnitude should not qualitatively differ from those of the standard scheme. 
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A Integrals of the strong interaction constraint 
diagrams 

This appendix is devoted to the calculation of the constraint diagrams entering in the 
strong interaction corrections to the potentials (Sec. 4 and Fig. 1). 

The incoming particles are designated by 1 and 2, the outgoing ones by 3 and 4, the 
loop particles by 5 and 6 . They represent either 7 r + 7 r - or 7 r° 7 r°. The relative momentum 
of the ingoing particles is designated by p, that of the outgoing ones by p'. The loop 


particles have momenta pi — k and P 2 + k. The constraint ( |2.1|) implies that pa, = Pl/2 
(i — 1 ,. .., 6 ) and k L = 0 inside the loop. 

The various momentum transfers have the following values: 

hi = (Pi - Pbf = k T2 , u 6 i = {pi - p 6 ) 2 = (2 p T - k T ) 2 , 

t 35 = {p 5 - P 3 f = ( P’ T ~P T + k T ) 2 , n 45 = (p B - P4) 2 = 0P' T + P T - k T ) 2 . 

(A.l) 


We also use notations (|2.9|) with b 2 = —p T2 and b' 2 = —p ,T2 . We choose for the evaluation 
of the integrals values of s below 4,m 2 , m n designating here the mass of the loop pion, so 
that 6 q is negative. The analytic continuation of \J—b^ to positive values of b g is done with 
the replacement \J—b^ —> — We use for the integrations dimensional regularization. 
We introduce the compact notation: 

(a%b) = j J^y a 9o(^J-,P T ~ k T )b, ( A - 2 ) 

where a and b represent either one of the momenta t and/or u defined above or the factor 
1. We find: 

d 3 k T 1 1 

47T 


So) = -/ 


(27t) 3 + ( P T — k T ) 2 + ie 

(dot si) = (gouei) = ]~(b 2 + 


-bo, 


(A.3) 


47T 

1 


(t 3 5go) — (^ 45 < 7 o ) — (b 12 + b^sj—bo, 


47T 

(t 33 got 3 l) = (w45i?oW6l) = — 


(t 3 5go u 6l) ~ ( u 45got 3 l) — 


47T 


§ + bv+ 5 bKb 2 + v 2 >+( bir 
\ubl + bV + ^(b 2 + bp + (blf 


bo, 

y/-%- 


(A.4) 
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Equations (|A.4| ) show that the integrals involving the variables t and/or u yield, with 
go.+~, quantities that are of order 0(a 3 ) and higher and thus are negligible. The similar 
terms considered with ( 70,00 yield contributions of order 0(a 2 \[^^), or 0((-^ 2L ) 3 / 2 ) and 
higher, which also are imaginary. Apart from being small, these terms will cancel the 
imaginary part of the amplitude present in the expression of the potential. The leading 
0(a) or 0(yj terms of the constraint diagrams are simply obtained by neglecting 
from the start the variables t and u in the scattering amplitudes, which then factorize 
outside the integrals and the latter reduce to the integral over the constraint propagator 
£0 [Eg. flOP]. 


B Two-loop diagrams 

In this appendix we generalize the validity of the lowest-order formula established in Sec. 
4 [Eq. (fi~7| )l to the two-loop level. 

The main ingredient of the proof is already met at the one-loop level and can be 
summarized as follows. In the vicinity of the two-pion threshold, the strong interaction 
unitarity one-loop amplitude (with the pion mass shift included in) can be decomposed 
into an analytic function Ad/ 1 / of the real variables 6 q ( s ), t and u and a non-analytic 
part, Ad// n (essentially proportional to \J— fejj). [The analyticity of Ad A in t and u near 
threshold is due to the absence of infra-red singularities in the strong interaction amplitude 
with massive pions.] The constraint diagram amplitude, Ad 0 *- 1 ’ 0,1 ) (the notation being 
explained before Eq. ( |5 ■ 8|) ), has the property of cancelling the nonanalytic part, A i^ n _, 
of the one-loop amplitude (sum of the two diagrams of Fig. 1): 

Ad (1) + M CW) = Ad/ 1 /. + M^ln. + -M c(1,0,1) = Ad/ 1 /.- (B.l) 

The deviation of the analytic piece, Ad//), from the pionium energy to the 7 r + 7 r - threshold 
is of order 0(a 2 ) and hence Ad A can be replaced by its value at threshold: 

M^l(b 2 0 ,t,u) = Ad//.(s = 4 ml+,t = 0,u = 0) + 0(a 2 ). (B. 2 ) 

Two cases must be distinguished, depending on whether the loop is made of 7 r + 7 r“ or 
of 7 r° 7 r°. In the first case ( 7 r + 7 r _ -loop), the nonanalytic piece, Ad/'/n > an d the constraint 
amplitude, Ad^ 1,0,1 ), are real and separately vanish at the 7 r + 7 r“ threshold and therefore 
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the value of M^n. threshold coincides with that of the one-loop amplitude A4 (1 ) at 
threshold, which is real. In the second case (7r°7r°-loop), and 7Vd c ^ 1,0,1 ' ) are imaginary 

and do not vanish at the 7r + 7r _ threshold; the value of Athen coincides with the real 
part of the one-loop amplitude Af A at threshold. Therefore, in all cases we have: 

M (1) + M C(1 ’0’ 1) = 7ZeM (1) (s = 4 m 2 +,t = 0 ,u = 0). (B.3) 

This property can be used for the analysis of two-loop diagrams. It is clear that not 
all two-loop diagrams necessitate a detailed study; those not having singularities in the 
s-channel are analytic in the vicinity of the two-pion threshold and real; since in this 
case their deviation from the threshold value is of order 0(a 2 ), they can immediately be 
replaced by their value at the 7r + 7r“ threshold. 

The typical diagrams involved in this analysis are presented in Fig. 4. 

The sum of the four diagrams of Fig. 4a (where the loop pions are the same in the 
four diagrams) is free of singularities in the s-channel and is represented in the vicinity 
of the two-pion threshold by a real analytic function. The value of the latter at the 
pionium energy differs from its value at the 7r + 7r“ threshold by an 0(a 2 ) term; hence 
it can be replaced by its value at the two-pion threshold. When the loops contain at 
least one pair of 7r + 7r - , this value coincides with the real part of the two-loop amplitude 
(first diagram) at threshold. If the two loops correspond to neutral pions, then this value 
differs from the real part of the amplitude at threshold by a factor that is proportional to 
the product of the imaginary parts of each loop (actually cancelled by the last constraint 
diagram) and hence to (Am T ) 2 . This factor contributes with a relative order of magnitude 
of 10 -4 and can be neglected. Therefore, the sum of the four diagrams at the pionium 
energy is equivalent to the real part of the two-loop amplitude (first diagram) at the 7r + 7r - 
threshold. 

In diagrams of Fig. 4b, the tadpole diagram factorizes at the vertex and does not 
interfere with the loop diagram calculations. For the latter, one has the same results as 
in the one-loop case (Fig. 1). 

In the first diagram of Fig. 4c, the internal propagator is modified with the inclusion 
of the self-energy correction (which must also be done in the other propagator). This 
feature does not qualitatively change the singularities of the diagram with respect to the 
one-loop case, for these depend essentially on the mass-shell condition of the internal loops. 
A similar self-energy inclusion in the constraint propagator [Eq. Q2.7D 1 (second diagram) 
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ensures the cancellation of the singularities of the first diagram and the reasonings of 
the one-loop case can be repeated. (The mass shift coming from the above self-energy 
correction can be incorporated in the mass term used in the constraint propagator go. 
Wave function renormalization constants should not influence physical results.) 

In the diagrams of Fig. 4d, the constraint diagram, corresponding to the first loop, 
cancels the singularities of the first loop of the first diagram. The second loop, which 
is on the right of each diagram, is free of singularities in the s-channel and is also free 
of infra-red singularities (no massless particles). Hence, the sum of the two diagrams is 
analytic in the vicinity of the two-pion threshold and the previous results are found again. 

There are also diagrams with three internal pion propagators, in which three external 
pions join each other at one vertex and a single external pion is attached at the other 
vertex. Such diagrams have singularities for pf > (3m f ) 2 (i = 1,... ,4). For pf s in the 
vicinity of the mass shell, as is the case in the present problem, these diagrams are free 
of singularities and are analytic in s, f, u. 

The remaining two-loop diagrams do not have singularities in the s-channel and the 
above diagrams exhaust the cases where constraint diagrams occur. The result found in 
Sec. 4 [Eq. (fT7|)] can therefore be generalized to the two-loop level. 

At the three-loop level, a qualitative change appears with the occurrence of inelas¬ 
ticities through the four-pion intermediate states. The constraint propagator (p.7|) is no 
longer sufficient by itself to cancel the singularities of the scattering amplitude. In this 
case, new pieces should be added to it to take into account the inelasticity effects. 


C Cancellation of linear divergences 


When the linear divergence of the continuous spectrum is kept, (Ar) conti [Eq. 
receives the contribution: 


® 3 )] 


(Ar) cont . — — 2r 0 H_i— ! _|— 


r°° dk 
L 27r 2 


+ 


(C.l) 


where the dots stand for the contributions other than that of the linear divergence. 

On the other hand, the integration of the constraint propagator g 0 [Eq. (|2.7|) 1 also 
produces a linear divergence. Equation ( A.3 ) now becomes: 


r d 3 k T 1 

J (27r) 3 &o + ( P T — k T ) 2 + ie 


dk 

2 ^ 



(C.2) 
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The corresponding linear divergence shows np in the contributions of constraint diagrams 
that enter in the calculation of the potentials. We first consider the constraint diagram 

^+-,+-<70.+-^+-, oo (hr lowest order) contributing to _,oo in Eq. ( |3.8|) . The linear 

divergence coming from go,+- [Eq. (|C.2|) 1 gives an additional contribution to E + _ i0 o : 

r°° dk 


14 -,oo — E + _ i+ _f / + _ 


00 


2 vr 2 


(C.3) 


This quantity, when inserted in the last term of the wave equation ( |3.13|) , provides, in 
first-order of perturbation theory, the contribution to the decay width: 

r°° dk 


PC 

(AT) c(2 ’°’ 2) = 2T 0 V + - +_ / 

Jo 


2t r 2 ’ 


(C.4) 


which cancels the linearly divergent contribution (|C.1|) . In the last term of Eq. ( |3.13|) 

there are also quadratic terms in f dk/(2ir 2 ), which however are proportional to V 2 , ,_. 

These must be associated with similar divergences arising from third-order perturbation 
theory. The constraint diagrams of Fig. 2 considered in Sec. 5 do not contain linear 
divergences and hence the net result of Eq. (|5.15|) remains unchanged. 

The other linear divergences present in the constraint diagram contributions must sim¬ 
ilarly be associated with other contributions having the same types of divergences. Thus, 
the linear divergence of the term V’ + _ j+ _< 7 o r i— V’+_ i 4 ._ which enters in the calculation of 
V+_ )+ _ [Eq. ( p.8|) j must be associated with the second-order perturbation theory contri¬ 


bution of V+_ !+ _ of the wave equation (|3.13[) which concerns the shift in the real part of 
the energy. The linear divergence of the term |E + _ j0 o^o,oof / oo,+^, contributing to 
cancels a similar divergence coming from the last term of the wave equation ( |3.13| ). The 
linear divergence of the term IVf^oof/o^oEoo.oo contributing to the calculation of I/ + _ 00 
must be associated with the contribution of the last term of Eq. ( |3.11|) (and which was 
neglected in the text because of the smallness of its finite part; cf. the paragraph preceding 
Eq. Q). 

In summary, the contributions of linear divergences, when isolated and grouped in 
terms with definite structeres, disappear by mutual cancellations from physical quantities. 
In this respect, the use of dimensional regularization automatically takes into account from 
the start the above cancellations. 
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D Integrals of the radiative corrections 


We present in this appendix the definitions and expressions of the integrals involved in 
the electromagnetic radiative corrections. For the evaluation of the divergent integrals we 
use dimensional regularization and designate by d and /j the space-time dimension and 
the regularization mass, respectively. 

The integrals involved in the self-energy part are the following: 
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fi?’(p) 
U 2 (p) 


= P 

= P 
i 


4 -d 


d d k 
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where we have defined: 
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One finds for the self-energy the expression: 


p rv 
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The integrals involved in the vertex function are the following: 
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Their expressions are: 
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One finds for the vertex function the expressions (at leading orders): 
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Pi A 2 / 

Among the integrals involved in the calculation of the constraint diagrams, F^ 0 and 
R C were given in Sec. 5 [Eqs. (|A9|) - (|5.11|) ]. The integral F^ c is: 
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E The scattering amplitude in the unequal mass case 


In this appendix we present the details of the calculation of the off-mass shell vr + (pi) 7 r _ (p 2 ) 
~ i> 7 r °(P 3 ) 7 r °(P 4 ) scattering amplitude at t — u — 0. We need for this evaluation the 
expression of the function R [Eq. ( [D . 1 0| )] in the unequal mass case. DeEning: 
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4 -d 
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with: 
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(E.3) 

In Eq. ( |E.2| ) is equal to £ [Eq. (p. 6 | )] with m\ replacing m w .. 

For small negative values of s, denoted now by t, R " behaves to first order in Am 2 = 
(m 2 — m|) as: 


R 12 (t) = - + 


1 . . 2 t. 


:(Am“ + -)- 
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The loop contributions with 7 r + 7 r _ and 7 r° 7 r° propagators, respectively, are: 


iF o M L P = -( s-2mB 0 ) 

= i(s-2A-B 0 ) 
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(E.5) 

(E.6) 


where T + and T° are dehned in Eq. ( |D . 1|) , with the 7r + and 7T° masses, respectively. 

The sum of the two contributions with a 7 r + 7 r° loop is, to first order in Am 2 (regular 
terms in t and u are neglected): 

lF o M ioo P = (Pi + P 2 + - 4mS 0 )T + + (pj + pj + 2m 2 0 - 4 mB 0 )T° 
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The tadpole diagrams at the four-pion vertex originate from the following part of the 
chiral effective lagrangian: 
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and yield the contribution: 


iF*M tad . = -[3sT° + AsT + -(s-pl-pl~2ml + )T + 
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The contribution to the four-point vertex function from the pseudoscalar densities 
(16.17| ) and (|7.5|) , concerning its mass shift part, is: 
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The 0(e 2 p 2 ) lagrangian (|6.16|) provides the contribution: 
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Figures 


Fig. 1. A loop diagram and its constraint diagram, denoted by a cross, in the s-charmel. 
The sum of the two diagrams is free of singularities in the .s-charmel. 

Fig. 2. Constraint diagrams the dominant part of which cancel the ultra-violet divergence 
of the continuum contribution in second-order perturbation theory. Wavy lines are 
photons. 

Fig. 3. One-photon exchange diagram in the presence of the strong coupling (a) and its 
constraint diagram (b). 

Fig. 4. Two-loop diagrams having constraint diagram counterparts. 
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